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In the context of a conformal Supergravity (SUGRA) model in the Einstein frame, in which 
the (next to) minimal supersymmetric Standard model can embedded naturally to produce chaotic 
inflation scenarios, we study properties of gravitino in the cases where it is stable or unstable. In the 
latter case, we demonstrate that for large dilaton scale factors there is an enhanced magnitude of the 
gravitino width, when it decays to neutralino dark matter, as compared with the Standard SUGRA 
case. In this context, we discuss the associated consequences as far as Big Bang Nucleosynthesis 
constraints and avoidance of gravitino overproduction are concerned. 



I. INTRODUCTION 



In m, [2[ simple classes of supergravity (SUGRA) mod- 
els describing superconformal coupling of matter to su- 
pergravity have been considered. The models con- 
tain non-minimal scalar/space-time-curvature couplings 
of the form $ R, where $ is a frame function, depend- 
ing in general on matter supermultiplets, including dila- 
tons. Such couplings have been argued to lead natu- 
rally to Higgs-inflation in both non-supersymmetric 
0] and supersymmetric theories 0-[ll ■ Scale- free globally 
supersymmetric theories, such as the Next to Minimal 
Supersymmetric Standard Model (NMSSM) [T^] with a 
scale-invariant superpotential, can be naturally embed- 
ded 0,01 into this class of theories, leading to new classes 
of chaotic inflationary scenarios Moreover, such mod- 
els have been considered in [ll| in connection with the 
possibility of dynamical breaking of supergravity the- 
ories, exploring further the conformal couplings of the 
gravitino four-fermion interactions. 

It is the point of this article to discuss the proper- 
ties of gravitino fields in such models, in particular in 
the context of the NMSSM. Specifically, we shall anal- 
yse decay processcs involving gravitinos and calculate the 
corresponding life time. Depending on the strength of 
the conformal couplings, the width can be suppressed 
or enhanced significantly. In the case of enhancement 
the rapid decay of the gravitino implies a resolution 
of the gravitino overproduction, avoiding the Big-Bang- 
Nucleosynthesis (BBN) constraints. 

The structure of the article is as foUows: in the next 
section |IT1 we describe the basic Lagrangian formalism 
and properties underlying the conformal supergravity 
models of In section Hlll we analyse the main decay 

processes involving gravitinos, calculate the associated 
widths (and lifetimes), and discuss how the latter are 



constrained by BBN. Conclusions and outlook are pre- 
sented in section ITVl 



II. LAGRANGIAN FORMALISM OF 
CONFORMAL SUPERGRAVITY MODELS 



The action of the conformal supergravity models of 
[l|, that we shall employ in our analysis below, in the 
Jordan frame reads: 



\<i> [Rie) - {^,R^] - i(5^1>)(Vi • lin + 



+£o + Cl/2 +Ci-V + Cm + Cnu^ + Cd + Ci{ , (1) 

where the curvature R{e) uses the torsionless connection 
ujf^"'''{e), with e^^ the vielbeins, and e the vielbein deter- 
minant, and the gravitino kinetic term is defined using 

i?^ = 7''"- {dp + ic^/''(e)7afc - ii^p75) • (2) 

Here Afi is the part of the auxiliary vector field containing 
only bosons, namely: 

- ii (a^2"a„/c - d^z^d^ic) - ^a^^Pa , (3) 

where A^^^ is the Yang-Mills gauge field, z" are (com- 
plex) scalar fields, /C(z, z) is the Kahler potential and Pa 
is a momentum map or Killing potential, which encodes 
the non-Abelian gauge transformations on the scalars 
and may also include Fayet-Iliopoulos terms. 

The notation Co, £1/2 and Ci denote respectively the 
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kinetic tcrms of spin 0,i,l fields in (P) ^: 

^1/2 = -k~9^pX^^x" + 



u B 



(Im /^s) F^^F^'' ^ + (a^ Im fAs) A^757''A^ 



(4) 



function $ in front of the four-gravitino terms. These 
imply that, depending on the value of $, assumed to be 
stabilised appropriately by rolling to the minimum of an 
appropriate dilaton potential (generated by e. g. by string 
loops, in oase one embeds such conformal SUGRA mod- 
els to string theory, or other ways of breaking the scale 
symmetry), the effective couphng of the four-gravitino 
interactions can be much larger than the gravitational 
coupling. Indeed, in the Einstein frame (denoted by a 
suffix E) , the graviton field (in the original Jordan frame, 
denoted by a sufRx J) is redefined by means of 



where fAsi^) is a holomorphic kinetic gauge matrix, 
La = daln{-^), La = La, L^p = daLp - and 
= ~l^9a0 + z^LaLp, with g^^ = dad^K. is the 
Kahler metric, with the notation da = gf^, da = gf^- 

In the notation of the covariant derivatives of the 
gauginos A"^ are defined as: 



{d. 



A 

(5) 



with Jab~^ the structure constants of the Non-Abehan 
gauge group. 

The fermion mass terms, including gravitino bare 
mass terms (if any) and the mixed terms £mix containing 
scalars and fermions, including factors of the frame func- 
tion, are given exphcitly in [ij, and again wiU not be of 
interest to us in this work. We shaU be exphcitly inter- 
ested in the penultimate of the terms on the right-hand 
side of Eq. ([T]), namely: 

U = i(Re/^s)V^^7"' [Fi -f f£) 7^A^ 

-i/ABar7"'i^ab^A^ - ici>L„r7'^-i?^7^, +h.C.} , 

(6) 



where 



r „h 



The explicit expression for the 4- fermion terms >C4f, 
which also contain a significant dependence on the frame 
function $ and its derivatives, will be presented be- 
low. Such four-fermion terms, in particular four-gravitino 
oneS; have been argued in to play an important role 
in some variants of the above class of conformal super- 
gravity models, which can characterise certain low-energy 
limits of superstring theories, in which the frame func- 
tion $ may be identified with the dilaton-axion complex 
superfield, ^ = ^ e"^*^. Such models can serve as proto- 
types in which the Deser-Zumino [l^ mechanism for dy- 
namical breaking of local supersymmetry (supergravity) 
scenario is realised explicitly. The important point to no- 
tice in these class of theories is the presence of the frame 



9 J =e ^9'^ 



«2 



9e 



(<I'|bo 



bosonic 



>o) , 

(8) 

so that the curvaturc term in the target-space supergrav- 
ity action has the canonical form, with coefficient the 
gravitational coupling = SttGn- with G n Newton's 
(four-dimensional) gravitational constant. 

In the Einstein frame, where we shall work for the pur- 
poscs of this work, the gravitational part of the effective 
conformal supergravity action, including the fermionic 
torsion induced four-gravitino terms reads 



^■^V T/S 



33 

64' 



s^^^ + p\x) + ^«p(:,)e-^ (^V^"^) + 



33 



(9) 



where R{e) denotes the cirvature term with respect to the 
torsion-free spin connection, denotes the canonically- 
normalised gravitino with Standard kinetic term as in 
= 1 supergravity, 

V'; = e'^V'M (10) 

while the . . . denote structures, including auxiliary fields, 
that are not of direct interest to us here. In writing ^ 
we have expanded the four-gravitino terms into detailed 
structures to exhibit explicitly the terms that generate 
masses, and we linearise the four-gravitino terms. The 
condensate of interest to us is the v. e. v. of the linearizing 
field p{x). 

The reader should notice that the coefficients of the 
gravitino-p interaction terms in (|9]) contain dilaton- 
dependent factors ~ e~"^, and are thus proportional, not 
to K, but to: 



K = Ke 



(11) 
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A one-loop analysis shows that the efFective potential for 
the condensate p field acquires a minimum at [TT| 

Pnun = {p) - ±0.726 (12) 

at which it vanishes. The gravitino mass term, thcn, in 
(jO]) then takes the fohowing form: 

- ™3/2'/''^r^''l/'',. = -^W3/2V''^l/''^, (13) 

with the dynamicaUy-generated gravitino mass of order 

2 408 

7713/2 = \/TTk- Vmin = 2.408 K-^e'f ^ —= Mp . 

v Stt 

(14) 

For large negative values of the v. e. v. ((p) < O the rc- 
suhing gravitino mass is mauch smaller than the Planck 
scale, and thus the efFective couphng R ([TT|) is much larger 
than the gravitational couphng Gn- 

This imphes that quantum gravitational corrections 
to the Minkowski space-time background, on which the 
above minimisation of the effective potential has been 
considered are not strong enough to destabilise (at least 
quickly) the gravitino condensate, unlike the case of Stan- 
dard N=l supergravity [isj . This prompted the authors 
of [ll| to consider large positive values of 

(e-'-^) = ^«'(^Ibosonic) » 1 , (15) 

and diseuss their relevance to the above-mentioned sce- 
nario of dynamical metastable breaking of local Siper- 
symmctry (SUSY) and gcneration of gravitino mass. 

It is the point of this article to examine the constraints 
implied by such an assumption on the dark seetor of the 
Universe in the case of the NMSSM embedded in this 
conformal supergravity framework. However, we shall 
not only restriet ourselves to the case of negative ex- 
pectation values of the dilaton, but we shall be more 
general and also consider the dark-matter phenomenol- 
ogy/cosmology of the case of positive dilaton v. e. v. In 
this latter case, dynamical breaking of local SUSY is not 
possible in view of the destabilising effects of the gravi- 
ton fluctuations, nevertheless one assumes conventional 
breaking of SUSY, e. g. through gluino condensation [l3 |. 
which is then communicating to the gravity seetor to re- 
sult in a non-trivial gravitino mass term TO3/2. 



III. DECAY PROCESSES INVOLVING 
GRAVITINOS AND ASSOCIATED 
CONSTRAINTS 

We first notice that, in the Einstein frame, one has to 
first normalise the kinetic terms of the scalars z" and 
X and A fermions, by appropriate redefinitions involving 
the frame function In particular the gauginos should 
be renormalised in the Einstein frame as 



in order to acquire canonical kinetic terms. On the other 
hand the gauge terms of the spin-1 part £1 (g]) are already 
in canonical form, in view of the conformal nature of 
the (Maxwell-like) kinetic terms for the Yang Mills fields. 
Here we consider the case /as= const = 1. 

Taking pU)) and ([T5)) into account, as well as the fact 
that the Yang-Mills gauge fields are not renormalised in 
the Einstein frame by the frame function e"^"^, we may 
write the first term of the right-hand side of © in the 
Einstein frame as: 

{Ld)E = i(Re/AB)<7"'' {f' t, + T'ti) . (17) 

This term is responsible for the gravitino-gaugino-gauge- 
boson interaction and it does not transformed going to 
the new frame. 

Tying to recap the transformations we have 

A' = e^^A. (18) 

The scalar and vector fields don't change. Using these 
we can calculate how the interactions that are relevant to 
the gravitino decays change due to the dilation presence. 
In particular we are interestcd in the ipxZ{'j) decays, 
where with x we denote the neutralino that is the LSP 
in our model. Doing so, we may then consider the terms 
Cd (H]) in order to compute the decay rate of the massive 
gravitino field G into , say, a neutralino x in the NMSSM 
and a Z gauge boson: 

G -^x + Z . 

The so affected gravitino decay rate, will in turn affect 
the Dark Matter relic density (assumed to be dominated 
by neutralinos in NMSSM) and this may imply stronger 
Big Bang Nucleosynthesis (BBN) constraints. It is there- 
fore important that detailed cosmological studies of such 
dilaton extended mSUGRA models are performed. 

In the usual mSUGRA gravitino satisfies the Rarita- 
Schwinger (RS) eqs. 

ii0-m)^^{x) = O (19) 
which result from the RS action [l^ 

- -^e'^''''"^M757.apV. - lmy,4'^.[r,ln^.^ ■ (20) 

The RS action including the dilaton effects can be written 
as 

>c' = -le^-'^i^'^^.jM - l^mi^'^b'^^W. (21) 

where -0^ — > G'^i^^i, and the relation between ttl'^i^ to 
TO3/2 is 



A' = e^^ A 



(16) 



TO3/2 = e'^m3/2 



(22) 
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so for ((/3) < 0; which is the physical case in several of 
the backgrounds discussed in order to allow for dynani- 
ical brcaking of local SUSY, the gravitino mass in the 
conformal supergravity scenario wiU be smaller than the 
corresponding one in the normal SUGRA. 

In the NMSSM the neutrahno field can be written 
as: 



where the spin average amphtude squared is 



l-^7P = T^ZI^ - K) (3^3/2 + K) (25) 



6Mf, m^i^ 



and the kinematical factor is defined as 



TVnS + N^2Wl + 7Vi3#(} + N^^Hl + N^^S , (23) 



J^{mo,mi,m2) = — j [("^0 ~ ("^1 + "^2)^ 



where N^j are the clemcnts of the 5x5 neutrahno diago- 
nahzing niatrix. If A^n dominates the sum J2i=i 5 ^ii — 
1, then the hghtest neutrahno is bino-hke. On the other 
hand, if is dominant then the neutrshno is singhno- 
hke. The latest data of the LHC experinients, indicat- 
ing a Higgs boson mass in the bahpark of 125 GeV p^ . 
combined with other experimental data from B-physics 
and dircct dark matter searches, seem to disfavor the 
singhno case [l7l |. Thus, in the foUowing it will be as- 



sumed that the hghtest neutrahno is mainly bino. In this 
case, the dominant two body decay channcls for gravitino 
are G — 7X and G — >■ Zx- Nevertheless, even in the 
singhno-hke case those channels, especiahy the 7%, are 
dominant, mainly due to the large available phase space. 
Therefore, our assumption that the light neutrahno is 
mainly bino is sufficicntly gcneric. 

We note in passing at this point that, as shown in (isj . 
a time-dependent (cosmological) dilaton (which can run 
with the cosmic time before BBN) can reduce consid- 
erably the neutrahno relic density, thereby increasing 
the cosmologically-allowcd available parameter space of 
SUSY even beyond the LHC reach. In this article we 
ignore such efFects, focusing on the gravitino intcractions 
exclusively, and assuming that the dilaton in our case 
has been stabilised to its vacuum expectation value at 
the scale of SUSY brcaking (or at least it is approxi- 
mately constant during a cosmological epoch). However, 
even for stabilised dilatons, in the Einstcin frame, the 
ncutralino-pair annihilation processes are affected by the 
conformal couplings, in particular are enhanced for large 
e'^ > \. Such enhancement may reduce the relic abun- 
dances already in the constant dilaton case. We postpone 
a comprehcnsive study of such cffects on gravitino decays 
and the neutrahno dark matter abundancc for a future 
study. 

Below we shall consider two cascs: one, in which the 
neutrahno is the stable Dark Matter candidate and the 
gravitino is heavier, thus unstablc, and the other, in 
which the gravitino is cosmologically stable and thus con- 
stitutes the dark matter candidate. We commence our 
discussion from the former case. 

In such a case, the formula for the decay width G — > 7% 
without the dilaton effects, reads as [191: 



X (mg - (mi - 7712)^)] 



2m1/2 



(26) 



On the other hand, taking into account the dilaton ef- 
fects, i. e. considering the corresponding process in the 
conformal SUGRA model in the Einstein frame, the cor- 
responding width becomes 



1 \M' P 
I = — J^(cm3/2, ' 0) 



IdiT cm 



3/2 



(27) 



where 



B'. 



X (6c4m2/2 + (c2 + l)m2) . 



(28) 



Above it was defined c — e'^. Noticc that putting c = 1 
we recover the result of Eq. ([24)) . The factor is rclated 
to the bino {B) and neutral wino (VFg ) components of the 
neutrahno, that is B^ = Nn cosOw + sin0w, where 
Ow is the electroweak mixing angle. 

For the channel G — ?> Zx in the Standard (dilaton-free) 
SUGRA the width reads as: 



Zx 



1 \Mz\ 

16 TT 7713/2 



-T{m3/2,m^,Mz) 



(29) 



where 



\Mz\^ 



BI 



6 A/2 7772^2 

127773/o771vM7 



37773/2 



^3/2 



Mi 



'3/2 



^1/2 {K - Mi: 



+ (7 



M|) 



(30) 



where Bz = —NusinOw + Ni2cos0w ■ In the case of 
dilaton the same width becomes 



1 \M,\ 

16 TT 7773/5 



-J"(77l3 / 2 , 777^,0), 



(24) 



1 iTW'zP 

= T7--7::^^^'''^V2^^x^Mz) , (31) 



16 TT C7(73/2 
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where 



Bl 1 



Ml 



6 (c 



+ (-6c* + 0^ + 1) m 



3/2 



1 mf 



the width for the conformal SUGRA case, where dilaton 
efFects are taken into account. An important observa- 
tion concerns the faet that for c < 1, which is the case 
where dynamical (metastable) breaking of local SUSY 
mass be possible, according to the argumcnts of [ll[ re- 
viewed above (c/, discussion leading to Eq. (fT4|)) the de- 
cay of gravitino to photons and ncutraUnos is kinemati- 
cally forbidden. This case will face important constraints 
from BBN which wiU be discusscd below. 



-6 



3 (-2c^ 



TO, 



(-2c4 + c2 + 1) 



1) 
1 



3/2 
2 



'■3/2' 



to: 



3/2' 



It is worth noticing that Tzx goes to F^^- 
AIz — >■ O, and the same holds for F'. 



(32) 
in the hniit 



1^ 




FIG 

G 

case. 
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1: The ratio 7? = F'/F for the channels 0—^X7 ^^nd 
XZ, as a. function of e''', for the neutrahno dark matter 



Wc start our numerical analysis discussing modcls 
where dark matter consists of neutrahnos and the grav- 
itino is unstable. In this case decays of the gravitinos to 
neutrahnos can be an important constraint, affecting sig- 
nificantly the BBN predictions. In Fig. [T]we present the 
ratio R = F'/F for the processes G — 7 x and G — > Z x 
as a function of e^. To make this figure we use the nu- 
merical values TO3/2 = 300 GeV, = 150 GeV, and for 
thesc the two body dccay widths for the dominant chan- 
nels (involving 7) are 1.4 x 10~^^ GeV and 3.3 x 10~^^ 
GeV for the 7 x Sind Z x, respectively. This yields a grav- 
itino lifetime ^ 4.7 x lO'' s. 

Wc first concentrate in the region of the Fig. [1] where 
c > 1, for which although dynamical generation of grav- 
itino mass and thus breaking of local SUSY may not oc- 
cur [nl, [1^ nevertheless, as mentioned above, one may 
assume a more or less conventional mechanism [l4j for 
SUSY breaking and the generation of a gravitino mass 
term TO3/2, (|2ip . (|22p . In this case, there is an enhance- 
ment of the ratio R = F'/F, where the prime denotes 




FIG. 2: The ratio R = F'/F for the channels x -)■ G7 and 
X ^ G Z, as & function of for the gravitino dark matter 
case. 

Conccrning the BBN constraints, we know that they 
becomc iraportant for lifetime of the unstable particle 
T > 10^ s [201 • Since t' = t/R and the gravitino lifetime 
without the dilaton effects is 0(10'') s, one observes that 
with i? > 10^ one avoids all the important BBN con- 
straints practically for any gravitino mass. Using Fig. [1] 
we understand that this happens for e*^ > 9 or 10. For 
smaller values of the BBN appears to become impor- 
tant, since the lifetime of the gravitino approaches its 
original value 10^ s as e"^ —s- 1, but its abundancc is still 
enhanced by the factor e'^ . On the other hand, values 
for e*^ < 1 cnhance dramatically the gravitino lifetime, 
leading to the exclusion of this range of the scale factor, 
being incompatible with the BBN constraints. We have 
checked numerically different values of gravitino mass up 
to TeV, but our results remain basically the same. That 
is, one finds that for e"^ in the range above 8 or 10, the 
BBN constraints can be avoided. 

In addition, we discuss also the complementary case 
where the gravitino is the LSP. In this case one stud- 
ics the reverse processes x ~^ G7 and x ^ G Z , from 
the point of view of the constraints induced by BBN. 
To compute the corresponding de cay widt hs we use the 
fact that the amplitudes squared [Al^^zP, both for the 
Standard and the dilaton cases, are the same as before, 
due to the assumed CPT invariance and unitarity. On 
the other hand, one has to interchangc and TO3/2 in 
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T and in the denominators of Eqs. ([24]), ([27]), (HH), and 
([gT]) . Doing so, we plot in Fig. [2]tlie ratio R = V' jV for 
these reverse processes. The numerical values we use are 
W3/2 = 150 GeV, 771^ = 300 GeV, exactly tlie reverse 
case of Fig. [1] With this choice, tlie neutralino lifetime 
is now ©(lO*") s. We thus see tliat the kinematicahy al- 
lowcd region happcns for values of c = e*^ < 2. This 
region depcnds on the choice neutralino and gravitino 
masses, but the attained values of R are not sensitive to 
this. One observes that in this case R < 10. Thus in this 
case the dilaton effects can not be used to relax the BBN 
constraints. 



IV. CONCLUSIONS AND OUTLOOK 



one has to assume morc-or-less conventional breaking of 
SUSY and its communication to the gravitational sec- 
tor. On the othcr hand, if the gravitino is cosmologi- 
cally stable, playing the role of dark matter, the BBN 
constraints are very restrictive in the full (kinematically- 
allowed) range of the scale (dilaton) factor e"^. 

We have not discussed in detail time-dependent dila- 
ton effects that involve running scale factors up to BBN, 
which are known to reduce significantly the dark matter 
abundances. Such effects, when combincd with the dila- 
ton effects on the gravitino decays considcred herc, may 
change significantly the cosmology and phenomenology 
of such conformal SUGRA models. We plan to return to 
these interesting issues in a forthcoming publication. 



In this article we have discussed the effects of sta- 
bilised dilatons on processes involving unstable particles, 
including gravitinos, that may affect BBN in conformal 
SUGRA models, incorporating the NMSSM in their spec- 
tra. We have found that, in the case where the gravitino 
is unstable, and the neutralino plays the role of dark mat- 
ter, there are regions of the scale (dilaton) factor e*^ > 8 
in which the BBN constraints can be avoided altogether. 
Moreover, in such regimes the neutralino dark matter 
abundances may be diluted thereby avoiding the cosmo- 
logical and particle physics constraints on SUSY mat- 
ter at current colliders, including LHC. Unfortunately, 
this case seems not to favour dynamical SUGRA break- 
ing (due to quantum-gravity instabilities) , and therefore 



Acknowledgements 

The work of N.E.M. is supported in part by STFC 
(UK) under the rcscarch grant ST/J002798/1 and by 
the London Ccntre for Terauniverse Studies (LCTS), us- 
ing funding from the European Research Council via the 
Advanced Investigator Grant 267352; this latter grant 
also supported visits by V.C.S. to the CERN TH Divi- 
sion, which he thanks for its hospitality and where this 
work was initiated. The work of V.C.S. is supported by 
Marie Curie International Reintegration grant SUSYDM- 
PHEN, MIRG-CT-2007-203189. 



[1] S. Ferrara, R. Kallosh, A. Linde, A. Marrani and 

A. Van Proeyen , Phys. Rev. D 82, 045003 (2010) 

[arXiv: 1004.07121 [hep-th]] . 
[2] S. Ferrara, R. Kallosh, A. Linde, A. Marrani and 

A. Van Proeyen , Phys. Rev. D 83, 025008 (2011) 

|arXiv: 1008.2942 [hep-th]] . 
[3] F. L. Bezrukov and M. Shaposhnikov, Pliys. Lett. B 659, 

703 (2008) [arXiv:0710.3755l [hep-th]]. 
[4] F. Bez rukov and M. Sh aposhnikov, JHEP 0907, 089 

(2009) [arXiv:0904. 15371 [hep-ph]] . 

[5] F. Bezrukov, A. Magnin, M. Shaposlmikov and 
S. Sibiryakov, JHEP 1101, 016 (2011) arXiv:100 '8H57| 
[hep-ph]]. 

[6] J. Garcia-Bellido, J. Rubio, M. Shaposhnikov and 

D. Zenhausern, larXiv:1107.2T63l [hep-ph] . 
[7] M. B. Einhorn and D. R. T. Jones, JHEP 1003, 026 

(2010) [arXiv:0912.2"7T8l [hep-ph]] . 

[8] L. AIvarez-Gaume, C. Gomez and R. Jimenez, Phys. 
Lett. B 690, 68 (2010) arXiv:100lT00^ [hep-th]]. 

[9] L. AIvarez-Gaume, C. Gomez and R. Jimenez, JCAP 
1103, 027 (2011) tarXiv:1101.4948, [hep-th]]. 
[10] C. Panagiotakopoulos and K. Tamvakis, Phys. Lett. B 
446, 224 (1999) [hep-ph/9809475] ; M . Maniatis, Int. 
J. Mod. Phys. A 25, 3505 (2010) |arXiv:0906.0"777l 
[hep-ph]]; For a review see: U. EUwanger, C. Hugo- 
nie and A. M. Teixeira, Phys. Rept. 496, 1 (2010) 



'arXiv:0910.1785' [hep-ph]] and references therein. 
[11] J. Ellis and N. E. Mavro matos, Phys. Rev. D 84, 085016 

(2011) iarXiv: 1108.08771 [hep-th]]. 

[12] S. Deser and B. Zumino, Phys. Rev. Lett. 38, 1433 
(1977). 

[13] L L. Buchbinder and S. D. Odintsov, Class. Quant. Grav. 
6, 1955 (1989). 

[14] see, for instance: H. P. Nilles, Phys. Rept. 110, 1 (1984) 

and references therein. 
[15] W. Rarita and J. Schwinger, Phys. Rev. 60, 61 (1941). 
[16] G. Aad ei al [ATLAS Colla boration] , Phys. Lett. B 716, 

1 (2012) [arXiv:1207.72T4l [hep-ex]]; S. Chatrchyan et 

al. [CMS Collaboration], Phys. Lett. B 716, 30 (2012) 
arXiv:1207.7235 [hep-ex]]. 
[17] K. Kowalska, S. Munir, L. Roszkowski, E. M. Sessolo, 

S. Trojanowski and Y. -L. S. Tsai. larXiv:1211. 16931 [hep- 

ph]. 

[18] A. B. Lahanas, N. E. Mavromatos and D. V. Nanopou- 
los, PMC Phys. A 1, 2 (2007) hep-ph/0608153 ; Phys. 
Lett. B 649, 83 (2007) hep-ph/0612152 ; A. B. Lah anas, 
Phys. Rev. D 83, 103523 (2011) ar Xiv:1102.4277 , [hep- 
ph]]; A. B. Lahanas an d V. C. Spanos, JHEP 1206, 089 

(2012) [arXiv:1201.260 1 [hep-ph]]. 

[19] J. R. EUis, K. A. Olive, Y. Santoso and V. C. Spanos, 

Phys. Lett. B 588, 7 (2004) hep-ph/0312262 . 
[20] M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 71, 



7 



083502 (2005) [astro-ph/0408426| ; R. H. Cyburt, J. EUis, 
B. D. Fields, F. Luo, K. A. Olive and V. C. Spanos, JCAP 



0910, 021 (2009) [arXiv:0907.5003l [astro-ph.CO]]. 



